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We accurately treat the efTect of the double and triple interactions of two-level systems (TLS) 
in glasses onto the energy delocalization due to the long-range interaction of TLS. Although this 
work qualitatively reproduces the estimates of our previous work we believe that it is important 
because in this paper the estimates are done with quantitative accuracy. This work can serve as the 
important step towards the development of quantitative theory of many-body delocalization due to 
the long-range interaction. Our study is compared with the recent work by Bodea et al [l[ claiming 
that the interaction of TLS triples leads to the energy delocalization. We cannot agree with the 
mentioned work because as we show the interaction of TLS triples was overestimated there and 
therefore the result for the TLS relaxation rate obtained there is invalid. This work serves as the 
extended complementary version of the one page comment submitted to Physical Review Letters. 
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The problem of Anderson localization in the presence 
of an interaction attracts a large attention since the con- 
cept of localization has been suggested in the classical 
work.Q Particularly, the long-range interaction l/i?" 
with a sufficiently small exponent a (a < 3) can lead to 
the delocalization at arbitrarily strong disordering. 0, Q 
The ensemble of two level systems (TLS) in amor- 
phous solids U is the particularly interesting system to 
study the effect of the long-range many-body interaction 
J(i?) ^ Uo/R^ on the localization, because the inter- 
action is extremely weak their, i. e. the dimension- 
less product of the density of two level system states Pq 
and their interaction constant Uq, is the universally small 
value around 1O~^.0 After the experimental demonstra- 
tion that the low temperature TLS relaxation rate shows 
linear temperature dependence 0,0, H instead of de- 
pendence due to the TLS-phonon interaction we propose 
the scenario, where this anomaly was explained by the 
TLS interaction stimulated delocalization leading to the 
irreversible relaxation. 

Recently the model of a TLS interaction stimulated 
relaxation has been considered by Bodea et al.[l| They 
obtained the estimate of the TLS relaxation rate differ- 
ent from Ref. 0] both qualitatively and quantitatively. 
The method used in Ref. [l[ was based on the analysis 
of resonant couplings developed in our previous work. Q 
The number of resonant couplings of single TLS with 
surrounding TLS (resonant pairs) and TLS pairs (reso- 



nant triples) were calculated. The resonant coupling was 
defined as in Refs. 0, namely two entities A and B 
arc in resonance, when their joint transition amplitude 
Vab induced by their interaction exceeds the difference 
between their energies Ea — Eb calculated ignoring their 
interaction. The criterion for energy delocalization was 
chosen using the condition of Ref. [9| that the given TLS 
have more than one resonant interactions. 

The number of resonant pairs W2 per the single TLS 
was found to be very small [H in full agreement with 
Ref. where this estimate was already made. Below we 
reproduce this estimate to save time of the reader. The 
Hamiltonian of two interacting TLS can be expressed as 

-A,5f-A,^-Ao.Sr-Ao,^-C/.,5^^. (1) 

Then we can develop analytical approach for double tran- 
sitions treating TLS interaction as a perturbation. To 
develop the perturbation theory with respect to the in- 
teraction we diagonalize the nonintcracting Hamiltonian 
introducing new spin operators 
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Then the Hamiltonian of pair can be rewritten as 
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The basis states of the Hamiltonian Hq can be repre- 
sented in terms of the projections of spins sf = ±s| = 
±1/2 =1 ± >. The transition ampHtude of the pair 
|H — >— >| — l->is defined by the as the corresponding 
matrix element of interacting Hamiltonian 

J i[/ ^21^. (4) 

The resonant pair is formed when 

\E,-Ej\<J,j. (5) 

The probability of resonance for the single pair can be 
estimated as g \ Jij \ , where g is the width of the distri- 
bution of a given TLS over energies. Then we need to 
sum all probabilities over all neighbouring TLS j. This 
sum can be expressed as the integral 

« 1 ^ 1 dr, I J,, I 
eiXr - r,)9{r, - ar) 
«PoC/olnf^) «1, (6) 

where the distance rj is taken in the range between their 
minimum value ar ~ (Uo/TY^^ and the resonant phonon 
wavelength At — c/T where c is the sound velocity. The 
distance cannot exceed the minimum distance because 
at shorter distances the interaction between TLS exceeds 
their thermal energy and therefore it is a very small prob- 
ability to find any of them in the excited state. The up- 
per constraint is because the interaction length is always 
restricted to the resonant wavelength. We also use the 
notation jph ^ Uo/{Xx^) corresponding to the TLS re- 
laxation time associated with their interaction constant 
as in [ij. 

However, in contrast with Q Bodea et al found that 
each TLS has the large number of triple resonant interac- 
tions with surrounding pairs of other TLS. According to 
[H this should lead to the delocalization of TLS energy 



and the irreversible relaxation of each TLS. On our opin- 
ion the estimate of the number of resonant triples W3 
made in Ref. [H exceeds the real value because the au- 
thors overestimated the triple TLS transition amplitude 
J assuming 

j' = iAjE.,fjf, + iA,/E,fjf,. (7) 



Here A and E = + Aq are asymmetry energy and 

energy of TLS, while Aq is its tunnelling amplitude. The 
expression Eq. ([7|) for the triple transition amplitude 
J leads to the result, which conflicts with the common 
sense, that if the TLS j is located arbitrarily far from 
TLS i and k {Rij « Rj^ 3> Rik) then the amplitude 
J does not depend on the long distance Rij and is de- 
fined only by the short distance Rik leading to the result 
J ~ UiJik ~ Uq/R^i^. On our opinion the triple TLS 
transition amplitude should vanish with vanishing the 
coupling between far separated entities. To prove that 
we perform analytical and numerical calculations of the 
triple transition amplitude and estimate the number of 
resonant triples that turns out to be much less then one. 

To calculate the triple transition amplitude we intro- 
duce the triple TLS pseudospin Hamiltonian 

3 3 ^3 

Htriple = - AiSl AoiSi - - ^ UijS^Sj. (8) 
i—l i—1 — 1 

The long-range interaction Uij « Uq /R^j essentially con- 
tributes to the delocalization for the far separated TLSs 
so we assume that the interaction is smaller than all TLS 
energies E = -^/A^ + Aq. This condition is satisfied in 
all previous studies for relevant TLS clusters. [J, Then 
we can develop analytical approach for triple transitions 
treating TLS interaction as a perturbation. To develop 
the perturbation theory with respect to the interaction 
we diagonalize the noninteracting Hamiltonian introduc- 
ing new spin operators similar to Eq. ([2]). The effective 
Hamiltonian reads 



Hi 



triplel 



Hq + U, Ho 



i=l 



U 



E, 



(9) 



Assume that we have resonant transition involving all 

three TLS between states s| = = —sf. = 1/2 < > 

^i = Sj — —s^ — —1/2 or cquivalcntly | + H — >< > 

h > . In resonance the transition energy Ei + Ej — Ek 

must be small so in our estimate we set Ei + Ej w Ek ■ 
The transition amplitude between states | + H — > and 



h > is defined by the TLS interaction U. In the 

first order in U one has Jijk + H — \U\ \- >= 0. 

The first term is equal zero because there is no three spin 
transition terms in the TLS Hamiltonian Eq. ^ because 
such terms must contain the product of three terms, 
while the maximum number of terms is two. In the 
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third order the transition ampUtude can be expressed as 
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where the sum is taken over all possible intermediate 
states of the Hamiltonian Hq Eq. ([9]). Taking all possible 
intermediate states we found the following final expres- 
sion 
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This result fully disagrees with the estimate Eq. ([7|) of 
Ref. When the distance between, say, TLS k and 

TLSs i and j becomes large the amplitude Eq. pT]) de- 



creases as Jijk oc i/Rii^ ^ ^/-^jk contrast with Eq, 



([7| where this amplitude is independent of the longest 
distance between TLS. 

We also verified our result for the triple TLS transition 
amplitude using exact numerical solution of the Hamilto- 
nian Eq. ([9]). The most straightforward definition of the 
transition amplitude is based on the repulsion of energy 
levels. If we have two nearly degenerate energy levels 
with energies Ei and E2 coupled by the perturbation 
matrix element V then the energy splitting of two levels 
can be represented as 
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Then the coupling matrix element can be defined as the 
minimum energy splitting for two energy levels. One can 
extract the transition matrix element V from the exact 
solution for the energy levels of the Hamiltonian Eq. ([9]). 
We minimize the splitting of energy levels with respect 
to the asymmetry energy of the third TLS A^. Gener- 
ally one can find four local energy minima corresponding 
to two pair resonances Ek « Ei and Ek « Ej and two 
triple resonances, which are realized if Ei + Ej w Ek and 
\Ei — Ej\ « Ek- We are not interested in pair resonances 
because their number has been already estimated and 
proved to be small. [l|,0] We study only one of the two pair 
resonances corresponding to the condition Ei + Ej Ek , 
which was assumed in our analytical estimate Eq. ([TT]) 
using the Scilab package, [l^. We have generated the pa- 
rameters of three TLSs using the random number gener- 
ator between and 1 which are = 0.544, Aj ~ 0.232, 
Aoj = 0.231, Aoj = 0.216 and Aofc = 0.883, take the 
displacement vector between TLSs i and k to be (0, 0, 2) 
and compute the transition amplitude as the function of 
the length of the vector connecting TLS i and j in the 
form Fjk = 77(200). The interaction has been set in the 



form Uij = 



1/4-; u,k = -1A4; u,k = 0.5/, 
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The 



graph in Fig. [T] shows the comparison of numerical ap- 
proach with our predictions Eq. pip and the predictions 



of Eq. dZD, Ref. d. It is clear that the approach of 
Ref. [l[ deviates from our estimate and exact calcula- 
tions by several orders of magnitude. The analysis of 
the second resonance leads to the same results. Thus we 
minimize the energy splitting between two energy levels 
of the Hamiltonian Eq. ^ assuming Ei + Ej w Ek- If 
we enumerate all eight energy levels in ascending order 
these would be fourth and fifth energy levels having the 
energies £^4^5 ~ i{Ei + Ej — Ek) ~ 0. The minimization 
leads to the resonant asymmetry energy A^*. 

One can estimate the number of triple resonances W3 
using the transition amplitude Eq. (|lip . This number 
can be obtained summing the number of all TLS triples 
involving the given TLS with the constraint that the en- 
ergy of TLS transition | Ei + Ej — Ek \ is less than the 
transition amplitude Jijk Eq. (jlip . This constraint can 
be approximately taking into account by replacing the 
integral over the third TLS energy with | Jijk |. Then 
the number of triple resonances per the single TLS i can 
be estimated as 
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This result is different both qualitatively and quantita- 
tively from the result of Ref. [ij . Remind that the con- 
tribution of such TLS results in the huge phase volume 
dimensionless factor which results in the large num- 
ber of triple resonances per one TLS in [l|. 

Thus according to our estimate, triples are not con- 
tribute significantly to the delocalization of TLS energy. 
This was already discussed in our previous work 0]. If 
however the more complicated clusters of four two level 
systems will be taken into consideration the number of 
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r~| numerical result (tliis work) 
^— analytical result (this work) 
^— estimate by Bodea et al 




Distance r., 



FIG. 1: Dependence of a triple tunnelling amplitude on the 
distance between TLS in different models. 
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Rmax ~- This four particle clusters made of resonant 
pairs of two-level systems are responsible for the energy 
delocalization as described in the previous work.Q 
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